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Superfluid phase transitions are discussed from a geometrical perspective as envisaged by Onsager. 
The approach focuses on vortex loops which close to the critical temperature form a fluctuating 
vortex tangle. As the transition is approached, vortex lines proliferate and thereby disorder the 
superfluid state, so that the system reverts to the normal state. It is shown in detail that loop 
proliferation can be described in exactly the same way as cluster percolation. Picturing vortex loops 
as worldlines of bosons, with one of the spatial coordinates interpreted as the time coordinate, a 
quantitative description of vortex loops can be given. Applying a rotation (to superfluids) or a 
magnetic field (to superconductors), which suppresses the formation of vortex loops and instead 
can lead to open vortex lines along the field direction, is shown to be equivalent to taking the 
nonrelativistic limit. The nonrelativistic theory is the one often used to study vortex lattice melting 
and to describe the resulting entangled vortex liquid. 



I. PREFACE 

As classical fluids and gasses, superfluids, Bose-Einstein condensates, and also superconductors support vorticesfl: 
extended line-like objects around which the fluid or gas executes a circular motion. In the classical world, vortices 
are better known as whirlpools (in fluids) and as whirlwinds (in the atmosphere) — or also as tornadoes, or typhoons, 
depending on which side of the Paciflc Ocean one lives. The velocity of the circular motion around a vortex line varies 
as the inverse of the distance to the rotation axis. The strong winds close to the center of a whirlwind — the so-called 
eye of a tornado — are, incidentally, responsible for most of the damages caused by typhoons. 

A striking difference with classical systems is that the circulation around these extended line objects in quantum 
systems is quantized. Given that a vortex usually is macroscopic in size and sets all the atoms or electrons of the system 
in motion, the quantization of circulation is truly remarkable and cannot in general be related to the quantization 
rules of quantum mechanics. Instead, as is well known, the quantization here is provided by topologyH. 

The presence of vortices in superfluids and BoseiJilinstein condensates is an emergent property, associated with 
the spontaneous breaking of global gauge symmetry^. In supercondutors, this symmetry becomes a local one, which 
strictly speaking cannot be spontaneously broken, but the following argument is not affected by it. The vacuum 
manifold, defining the system's minimum of energy, is a circle S^^^. parametrized by the phase of the spontaneously 
broken U(l) gauge symmetry. If a vortex is circled once, the phase moves around the circle S^j^j. representing the 
vacuum manifold. A vortex therefore provides a map from the circle in real space to S^j^j.. Such maps are 




FIG. 1: Drawing, introduced by the psychologist J. Jastrow (1900), that depicts both a duck and a rabbit. 
* Write-up of lectures presented at the Kevo Winter School, Kevo, Finland, 20-26 April 2002. 
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characterized by an integer winding number telling how often the map wraps around the circle 8^^^^. when S], in real 
space is circled once. The integer nature of the winding number explains the quantization of circulation around 
vortices in quantum systems. A vortex of a given winding number cannot be continuously deformed into one of a 
different winding number and hence is topologically stable. 

In these lecture notes, superfluid phase transitions are discussed from the perspective of quantized vortices in the 
spirit of Onsager's picture of the A transition in liquid ^Hco. The description envisaged by Onsager is one entirely in 
terms of these one-dimensional objects, with their geometrical properties such as fractal dimension and configurational 
entropy. In the absence of an applied rotation (for superfluids) or magnetic field (for superconductors) , vortices cannot 
terminate inside the system and must form closed loops. A superfluid phase transition is characterized in this picture 
by a fundamental change in the typical loop size. Whereas in the superfluid phase only finite loops are present, 



at the critical point infi* 
phenomena at criticalljol 



loops appear — similat-|to the sudden appearance of a percolating cluster in percolation 



In Feynman's wordsH "The superfluid is pierced through and through with vortex line. 
We are describing the disorder of Helium I." 

The similarity between vortex loop proliferation and cluster percolation is, as will be discussed, very close. It derives 
from a similarity in the cluster size and loop size distribution. Both have the same form containing two factors, one 
related to the entropy of a given cluster or loop configuration, and the other related to the Boltzmann weight assigned 
to the configuration. Close to criticality, both of these factors are parametrized by a single exponent. One specifies 
the algebraic behavior of the distribution at criticality, while the other describes how the Boltzmann factor tends to 
unity upon approaching the critical point. Clusters and loops that are otherwise exponentially suppressed can grow 
without bound when the Boltzmann factor becomes unity and, thus, gain configurational entropy without energy 
cost. 

The picture of the A transition in liquid "^He in terms of proliferating vortices, which disorder the|S¥stem and thereby 
restore the spontaneously broken U(l) symmetry, has been advanced since then by various authorsli3. Early numerical 
evidence for this picture based on Monte Carlo simulations of the three-dimensional XY-model was given IPtR efs. 
|l2|,|l3[ Analytic methods io describe the transition using vortex loops were further developed by Williamalil, and by 
Shenoy and coUaboratorscJ. More recent numerical work on the three-dimensional XY-modcl from the perspective of 
vortices, in particular their loop size distribution, can be found in Rcfs. |T6| , pT|JT^ , [l9| . 

The famous drawing (see Fig. |l|) introduced by the psychologist J. Jastrow serves as a metaphor for the approach 
followed in these notes. The drawing — a three-dimensional snow version of which was sculptured by students attending 
the Kevo Winter School to guard Vortexland — stands for superfiuid phase transitions, which can be described either 
in the conventional (duck) wayEI, or in the Onsager (rabbit) way pursued here. 

The next section discuses cluster percolation — the paradigm of geometrical phase transitions — in its pure form 
known as uncorrelated percolation (Sec. O) a nd i n its applied form known as c orre lated percolation to describe 

(Sec 



thermal phase transitions in spin models (Sec. II B) and the liquid-gas transition (ttlQ). Section III gives a parallel 



treatment of random walks, in its pure form (III A) and in its applied form to describe random (illB) and correlated 
vortex tangles featuring in sup erfluid phase transitions (IIIC) and also in superfluid turbulence and defect formation 
after a rapid phase transition ( III D ). The section is set up such as to highlight the close similarity between cluster 



describes vortex lines in deta il, starting with noninteracting vortices 
:VB|) and an external field (|lVC|) . Section discusses entangled vortex 



percolation and loop proliferation. Section 
( [ V A ) , followed by including interactions 
lines as they appear when the Abrikosov fiux lattice mel ts (VA) from the perspective of Feynman's cooperative 
exchan ge ri ng theory of Bose-Einstein condensation (V_B) and also discusses the order parameter describing this 
phase (|VC|). 



II. CLUSTERS 



In this section, the paradigm of geometrical phase transitions, viz. percolation, and its extension to describe thermal 
phase transitions are discussed. 



A. Uncorrelated Percolation 



To define (site) percolation, consider a lattice with N sites labeled by the index 1 < i < iV. Imagine visiting a site 
and generating a random number < < 1. If that number is smaller than a predefined number p, the occupation 
probability, occupy the site, otherwise leave it unoccupied. This algorithm can be summarized in pseudo computer 
code as: 



while i < N 



FIG. 2: Typical outputElJ of the site percolation algorithm implemented on a two-dimensional square lattice for p < pc (left 
panel) and p ~ Pc (right panel). Occupied sites are marked by a dark square, with those belonging to the percolating cluster 
having a lighter grayscale, while unoccupied sites are unmarked. 



generate < < 1 
if ri < p occupy site 
i + + 
end while 



After all lattice sites have been visited in this way, a random pattern of occupied and unoccupied sites emerges (see 
Fig. Next-neighboring occupied sites can be grouped together into clusters, whose properties as a function of p 
form the subject of percolation theory. Since a given site is occupied whether or not neighboring (or any other) sites 
are, this process is more precisely referred to as uncorrelated percolation. 

As the occupation probability p increases, more sites become occupied and the clusters grow. For high-enough 
values of p, a large cluster, spanning the lattice is frequently generated during updates (see the right panel of Fig. y). 
The appearance of such a percolating cluster is a fundamental property and not an artifact of working on a finite 
lattice. Even on an infinite lattice, a percolating cluster, which would then be of infinite size, appears. The only 
difference between a finite and infinite lattice is that the percolating cluster on an infinite lattice first appears at a 
specific value Pc, special to the dimensionality and symmetry of the lattice considered, whereas on a finite lattice there 
is always a finite probability to generate a percolating cluster when p < pc. Also, for p > pc, the percolating cluster 
is always generated on an infinite lattice, but not necessarily so on a finite one. The theory of finite-size scaling is 
needed to extract the-i'alue of pc from numerical simulations done on finite latticesEll. 

Percolation theorycj can be set up starting from the cluster size distribution £s{p), giving the number density of 
clusters of size s. Close to the percolation threshold, it takes the form 

4(p)(xs-"e-^^ c^{p,-p)^/^, (1) 

where the coefficient c vanishes with an exponent 1/cr when the percolation threshold Pc is approached from below. 
The cluster size distribution (|l|) contains two factors. The first one measures (as will become clear below when 
discussing random walks) the configurational entropy of clusters, while the exponential is similar to a Boltzmann 
factor which for finite c suppresses large clusters. The latter vanishes when the percolation threshold is approached 
from below, implying no longer a restriction on forming arbitrary large clusters and resulting in a proliferation of 
clusters. 

Although percolation lacks a conventional description in terms of an Hamiltonian, a partition function Z can 
nevertheless be defined through 



ln(Z)ocV^4, (2) 



with V denoting the volume of the system. The right hand counts the number of clusters of all sizes. As in statistical 
physics, various physically relevant quantities can be calculated from Z , such as the percolation strength P{p), which 
for p > Pc denotes the probability that a randomly chosen site belongs to the percolating cluster, and the average 
cluster size S{p). Specifically, 

oc oc/ 

showing that the percolation strength P is similar to the magnetization in spin models, while the average cluster 
size S is similar to the magnetic susceptibility. Near the percolation threshold, the quantities (^ show power-law 
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behavior, P{p) ^ {p — Pc)^^°^ and S{p) ~ \Pc— p\~^'^'"' i characterized by the critical exponents /3por and 7pcr, where the 
subscript "per" is to indicate that they pertain to percolation. A final critical exponent, ?7perj specifies the power-law 
behavior of the correlation function Gper at the percolation threshold. Physically, GperCx, x') denotes the probability 
that the sites with position vectors x and x' belong to the same cluster. It usually is only a function of the distance 
between the two sites and has &t p = Pc the algebraic behavior in d space dimensions 



Gpcr(:r:) - pr^- (4) 



The exponent tr, together with the so-called Fisher exponeninJ r, specifying the cluster size distribution ([l]), determine 
these critical exponents through scaling relations [see Eq. (|l^) below] . 

Since the percolating cluster, often also referred to as infinite cluster, includes a macroscopic fraction of the total 
number of sites, its mere size would dominate any sum. For this reason, summations over cluster sizes exclude 
the percolating, or infinite one. In practiceEil, when working on a finite lattice, simply the largest cluster present 
is excluded, whether or not is spans the lattice. The percolation strength P can nevertheless be extracted from 
considering only finite clusters in Eq. ( p^ ) because of the constraint 

p^p{p) + Y,-^Up). (5) 

stating that an occupied site either belongs to the infinite cluster or to a finite one. This constraint is special to 
percolation. 

Another important quantity characterizing percolation is the radius of gyration Rs , 

i?^ = i^(x,-x)^ (6) 

which gives a measure of the spatial extent of a cluster of size s centered at x = (1/s) X]i=i '^i- For large enough 
clusters, Rs scales with s as 

Rs - s^'^, (7) 

defining the Hausdorff, or fractal dimension D. Close to the percolation threshold, D can be related to thCjCorrelation 
length exponent specifying how the correlation length ^ diverges, £,{p) ^ \pc — p\~" . The relation readsE3 

Also the Fisher exponent r can be easily related to the fractal dimension, giving 

^ = ;^ + 1, (9) 

where d denotes the number of space dimensions. 

Using scaling laws, one can express all the pjtical exponents in terms of just the two exponents a and t, specifying 
the cluster size distribution (|^). For examplec3, 

T — 2 3 — T 

/3pcr = , 7pcr = , ?7por = 2 + d - 2D, (10) 

cr cr 

with a increasing from 36/91 « 0.40 on a two-dimensional square lattice to 1/2 on a six-dimensional square lattice, 
while at the same time, the fractal dimension increases from Z? = 91/48 « 1.90 to 4. Six is the upper critical dimension 
of uncorrelated percolation, beyond which the exponents remain fixed at their mean-field values tr = 1/2 and D = 4E3. 

The discussion so far concentrated on clusters of occupied sites, which for p < pc make up the minority phase, while 
the unoccupied sites form the majority phase. At the percolation threshold this changes and the roles of the minority 
and majority phases are interchanged. In other words, when approaching the percolation threshold from above, one 
would naturally concentrate on the disfavored unoccupied, rather than on the abundant occupied sites. From that 
side, the transition corresponds to a proliferation of clusters of unoccupied sites. In a sense, the two descriptions are 
dual to each other with the minority and majority phases interchanged. 
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B. Correlated Percolation 



To see how percolation theory can be apphed_±o thermal phase transitions taking place in spin models consider the 
Ising model at the absolute zero of temperatureEj. Here, in the absence of fluctuations, all the spins, which in principle 
can take the values 5'^ = ±1, point in one direction. When the temperature T increases, some of the spins will reverse 
their direction. The resulting spin configuration at a given temperature is easily mapped onto a percolation problem 
by identifying the sites with reversed spins with occupied sites. The question is now if the critical behavior of the Ising 
model can be recovered, focusing exclusively on these clusters instead of the magnetic properties. In other words, are 
the correct values for the critical temperature and the critical exponents obtained when measuring quantities like the 
percolation strength and the average cluster size discussed above? 

It turns out that, in general, clusters of sites with reversed spins proliferate already before the thermal phase tran- 
sition is reached. (Only on a two-dimensional square lattice, the temperature where a percolating cluster first appears 
coincides with the critical temperature.) Moreover, the critical exponents thus obtained are those of uncorrelated 
percolation and not of the Ising model (even on a two-dimensional square lattice). This means that the clusters 
constructed in this naive way lack-essential information about the Ising model and need to be modified. 

To construct a modified clustered, take two nearest-neighboring sites i, j of a naive cluster with its spins reversed, 
and include the pair into the modified cluster only with the bond probability 

= l-exp(-2/3J), (11) 

where J is the spin-spin coupling of the Ising model and (3 — with /cb Boltzmann's constant. The factor 2 J 

appearing here corresponds to the increase in energy when one of the two spins involved in the bond is fiipped. The 
upshot of this modification is that the resulting bond clusters are in general smaller than the naive ones and also more 
loosely connected. The modified clusters turn out to proliferate right at the critical temperature and display critical 
behavior identical to that of the Ising model — as should be accordiag to an exact map of the partition function of the 
Ising model onto correlated percolation by Fortuin and KasteleyrO. 

The above scenario was recently extended to three-dimensional 0{n) spin models, with n = 2,3,4. Instead of 
taking only the two values Si = ±1, which can be thought of as comprising a zero-dimensional sphere, the spins of 
the 0(n) model can take any value on the (n— l)-dimensional sphere defined by Sf = 1. The 0(n) model with n — 2 
corresponds to the XY-model, while that with n = 3 corresponds to the Heisenberg model. An exact analytic map 
onto correlated percolation as existed for the Isins-^odel is not available for n > 1, so that numerical simulations 
are needed to establish the existence of such a mapEj. Two nearest-neighboring sites i and j with spin S,; and Sj are 
included in a cluster with bond probability 

p^^ = 1 _ e-2/3Jmax[0,(S,.r)(S,.r)]^ (^2) 

where r is a randomly chosen unit vector which can be used to cut the unit hypersphere = 1 in two halves (see 
Fig. H). When the spins and Sj take values in opposite halves of the hypersphere, their projections onto r have 
opposite sign. Such a pair is never included in a cluster since pij = then. When the spins take values in the same 
halve, the pair is included in a cluster with the bond probability determined by the projection Si • r of their spin 
onto the random vector. The clusters thus constructed turn out to proliferate at the critical temperature where the 
thermal phase transition takes place. Moreover, the critical exponents obtained by measuring the percolation strength 
and average size of these clusters are the-aame as those obtained in the more traditional manner by focusing on the 
magnetic properties of these spin modelsBs. 




FIG. 3: A 2-sphere denoting the possible values a spin of the 0(3), or Heisenberg model can take. A randomly chosen vector r 
(not shown in the figure) can be used to cut the sphere in two halves, depending on the sign of the projections of a spin onto 
that vector. 
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FIG. 4: Sketch of part of the hquid-gas phase diagram, showing the coexistence curve and the critical point at which this curve 
terminates. 



These examples illustrate that thermal phase transitions taking place in a wide range of statistical physics models 
can be reformulated in terms of percolation theory. Each of these models is thus characterized by a partition function 
of the form with a cluster size distribution given by Eq. (|]). The exponents a and r characterizing take 
values specific for the universality class to which the model belong. Given the values of a and r, which in general 
differ from those of uncorrelated percolation, the critical exponents follow from scaling relations, as in Eqs. (p|)-(]lO|). 



C. Liquid-Gas Transition 



Next, the condeusation of a classical gas, the liquid-gas transition, is considered from the perspective put forward 
by Mayer et aiEMn and shown to be closely related to correlated percolation. 

The relevant part of a typical phase diagram is sketched in Fig. ^. Below a critical temperature Tc, the gas phase 
is separated from the liquid phase by a discontinuous phase transition. At the phase boundary, the two phases have 
the same free energy and despite different particle number densities coexist, separated by an interface. The energy 
of the mixed state is higher than that of the pure states because the formation of the surface requires energy — the 
surface tension. The coexistence curve ends at a critical point where the surface tension vanishes and the transition 
becomes continuous with Ising model critical exponents. The order parameter, distinguishing the liquid from the gas 
phase, is the difference in particle number density. i— . 

In Mayer's theory, which was the first description of a phase transition beyond mean-field theorycJ, the full interac- 
tion between molecules is approximated by a pair interaction. That is, three-particle and higher interactions, which 
are irrelevant when describing critical behavior, are ignored. Tlie .grand-canonical partition function Z{P,a), with 
a = —(3fj,, can then be written as an Ursell-Mayer cluster seriesEZru 

lr,iZ)=Yj2^e-"\ (13) 



where = h\j2'n:(3 /ra^ with m the mass-.of a molecule, is the de Broglie thermal wavelength, and the Ws, with 
&i = 1, denote the Mayer cluster integral£a. Being determined by the pair potential, these integrals can in principle 
be calculated, but, apart from small clusters, are in practice too complicated to evaluate for realistic choices of the 
potential. A (classical) ideal gas corresponds to setting 6; = for Z > 1. With the abbreviation 

^i = ^e--\ (14) 

the partition function (^3|) can be written in a form similar to Eq. (H) for percolation, 

HZ)^yY,Iu (15) 



while the expression for the number N of molecules, at least in the gas phase, becomes 

d_ 

da 



7V--Ain(Z)-V^Z£z. (16) 
I 
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It is tempting to interpret as the number density of clusters containing I molecules. However, the Mayer cluster 
integrals bi and therefore can be negative because the intermolecular interaction u{r) is attractive for large separation 
r, behaving asymptotically as u{r) ^ — In the gas phase, where large clustejca are exponentially suppressed so 
that only small clusters exist, the above interpretation is not obstructed too muchtj. 

As for percolation, a possible infinite cluster, containing a macroscopic fraction of the total number N of molecules 
present, is excluded in summations over cluster sizes. (In numerical simulations, simply the largest cluster is best 
excluded, whether or not it spans the lattice.) 

Condensation, achieved, say, by reducing the volume at fixed N , corresponds in this description to the formation 
of an infinite cluster. The sum / ti in Eq. ( p^ ) denotes the number density of molecules in the gas phase, which is 
represented by the finite clusters. When the coexistence curve is reached at the volume Vco+ , the sum saturates at the 
value N/Yco+- Upon reducing the volume further, an infinite cluster, which is not accounted for in the summation 
over cluster sizes, must form to accommodate the difference between the true density N/Y and that contained in the 
gas phase, N/Yco+- For this to be possible, the exponential suppression of large clusters needs to vanish. Since the 
chemical potential remains negative at the coexistence curve, the Mayer cluster integrals must yield an /-dependent 
exponential factor to cancel the factor exp(— a/) in Eq. ( [l^ at that curve. That is, 

^cx/-e"-', (17) 

close to the coexistence curve, with the first factor measuring the configurational entropy of clustersEl. It then follows 
that 



r^e-("-"-)', (18) 



so that as a approaches Oco from above, ii displays algebraic behavior, as required. This limit also marks tha.radius 
of convergence of the sum (p^: for a smaller than aco, the sum would diverge, i.e., at the coexistence curvetil 

A|e" = lim(/60^^'- (19) 

The discontinuous liquid-gas transition fundamentally differs from the continuous phase transitions discussed so 
far. In the earlier discussion, the appearance of an infinite cluster signaled the (continuous) transition of the entire 
system from one phase into the other. That is, the infinite cluster together with the smaller clusters still present (see 
the right panel of Fig. |^) represent the new phase. In the liquid-gas transition, the situation differs in that the infinite 
cluster denotes the liquid phase, while the smaller clusters denote the gas phase. In other words, the right panel of 
Fig. ^ now represents a snapshot of a two-dimensional section of the system (in microgravity conditions, say) with 
both the liquid phase (percolating cluster) and the gas phase (smaller clusters) coexisting. 

This interpretation is suppprted by the observation that the pressure associated with the finite clusters satisfies the 
Clausius-Clapeyron equationEll. To first order, the pressure P = ln(Z)//3V in the gas phase can be approximated by 
that of an ideal gas, corresponding to the first term in the Ursell-Mayer expansion (p3|). When differentiated with 
respect to T, this gives 

dln(P) 1 ,. I dln(b0 



where use is made of Eq. (19). To interpret the last term, the energy E = — dln(Z)/d/3 of the gas phase is determined. 



where a is to be considered a /3-independent variable, leading to the result 

E = yY,^M, (21) 

with e; the energy per molecule of a cluster of size /, 

..^IfcBT + fcBT^l^. (22) 

Since d&;/dT is in general negative, e; is smaller than the average kinetic energy per molecule, ^k^T. Physically, the 
last term in Eq. ( p^ ) denotes minus the (asterage) binding energy vi per molecule. In terms of this quantity, Eq. ( [20| ) 



becomes the Clausius-Clapeyron equatiorO 



dln(P) _ kBT + v^ 
dT " kBT^ ' ^ ' 
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FIG. 5: Schematic representation of the radius of gyration _R„ of a closed random walk containing n steps. 

with the numerator at the right hand denoting the energy required to vaporize one molecule from the liquid, viz. the 
sum of the thermal energy /cbT and binding energy Voo- 

It should be kept in mind that the liquid-gas transition, being discontinuous, is special in that both the liquid and 
gas phase can coexist, each occupying separate parts of the total volume. Similar arguments as the ones developed 
for the liquid-gas transition are sometimes applied to an ideal Bose gas to argue that Bose-Einstein condensation is a 
discontinuous phase transition. A Bose-Einstein condensate is, however, formed in the zero-momentum state invading 
all of space and therefore intertwined with the non-condensed part. No interface separating the two parts exist. The 
situation is similar to that in the 0{n) spin models discussed above, where, at the critical temperature, the entire 
system transforms into a new phase that cannot be spatially resolved into two separates parts, each occupying its 
own volume. 

Recently, the three-dimensional threa-state Potts model in an an external field has been numerically investigated, 
using clusters as fundamental objecta23. The phase diagram of this model is very similar to that of the liquid-gas 
transition, having a line of discontinuous phase transitions at small fields that terminates in a point at some critical 
field value where the transition becomes continuous with Ising model exponents. 

III. RANDOM WALKS 

The discussion of phase transitions so far was modeled after percolation theory, with clusters specified by the cluster 
size distribution (|l|) playing the central role. Another paradigm of thermal phase transitions is random walks to which 
the discussion now turns. 

A. Brownian Random Walks 

The simplest random walk is the Brownian random walkH. On a d-dimensional square lattice, where each site 
has 2c? nearest neighbors, a Brownian random walker located at a given site i randomly chooses among its nearest 
neighbors (n.n.) mi G {±j} the next site to visit. Here, the d unit vectors spanning the square lattice are denoted 
by j. The random walker is not subject to any restriction and can, for example, return to a site it visited before. In 
pseudocode, the algorithm describing a Brownian random walk reads: 

while i < n 

generate € {±j} 

move to n.n. rm 

i + + 
end while 

A central concept in the theory of random walks is the probability Kn{x — > x') for a random walker to move from 
one site with position vector x to another x' in n steps. The probability is given by 



9 



where the denominator, counting paths with arbitrary endpoints (indicated by is included so that Kn is a genuine 
probabihty, taking values between and 1. To evaluate it, imagine the random walker taking a first step to one of 
its nearest neighbors ±j. It then has only n—1 steps still available to reach the final destination x', with probability 
Kn-i{x + aj x'). That is, Kn satisfies the recurrence relation: 

if„(x^x') = ^^i^„-i(x + aj^x'), (25) 
±j 

where the factor 2d denotes the number of nearest neighbors on a d-dimensional square lattice. When _R'„_i(x x') 
is subtracted from both sides of Eq. (E5|) , a difference equation is obtained 



if„(x ^ x') - i^„_i(x ^ x') = ^ ^ [i^„-i(x + aj ^ x') - if„_i(x ^ x')] , (26) 

±j 

with the obvious boundary condition ifo(x x') = 5x,x'- The continuum limit is taken by simultaneously letting 
a and n —> oo, such that the length L = na of the random walk satisfies La — nc? — > const. The difference 
equation then turns into the differential equation 

a„i^„(x x') = ^V2i^„(x ^ x'), (27) 



with the celebrated solution 

/2 



Kn{y^ x') = 2 



(28) 



where no? — > const. Each dimension the random walker is free to roam gives a contribution ^Jdl^/rma? to the 
prefactor. 

The radius of gyration i?„ of a random walk containing n steps is defined by 



1 " 

Rl^-Y.^^^-^f■ (29) 



n 

1=1 



As for clusters, it measures the spatial extent of the random walk containing n steps and centered at x = (1/'^) Xi 
(see Fig. I). 
It scales with the number of steps as 

Rn ~ n^/^, (30) 

where for Brownian random walks, the fractal dimension is _D = 2. 

When the random walk is taken as representing a flexible polymer or vortex, the energy to create such a line object 
needs to be included. Let 9 denote that energy per unit length, i.e., the line tension, then this is achieved by the 
modification 

i^„(x ^ x') ^ X„(x ^ x') e-''«"^ (31) 

Because of the energy cost, the Boltzmann factor exponentially suppresses long lines, provided the line tension is 
positive. The total probability for a random walker to move from one site x to another x' using an arbitrary number 
of steps defines the correlation function G(x, x'), which is assumed to depend only on the distance between the two 
sites: 

G(|x-x'|) = ^X„(x^x')e~'5»-. (32) 

n 

In the absence of external sources and boundaries, vortex lines cannot terminate inside the system and must form 
closed loops. The partition function describing a (grand canonical) ensemble of closed, noninteracting vortex lines is 
given by 

ln(Z) = VV -X„(x x) e-'5^"", (33) 
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where the notation i4r„(x x) indicates that the random walk representing the vortex hne starts and ends at the 
same point. The right hand of the equation is therefore x-independent, as required for a partition function. Since a 
cfosed foop of length na can be traversed starting from any of the n sites visited, the factor 1/n prevents overcounting. 
With the explicit solution (^8|), the partition fimction becomes 

ln(Z)=V^£„, (34) 

n 

where 

£„cxn-"e-''''"", r=^ + l. (35) 

Physically, denotes the loop size distribution, giving the number density of loops containing n steps. In this way, 
the partition function is represented in a form surprisingly similar to that in percolation theory Eq. (||) and for the 
liquid-gas transition Eq. (Esl), with the proviso that not clusters, but loops are considered. 



B. Random Vortex Tangles 



Given the close similarity in describing uncorrelated clusters and random vortex loops, the question arises whether 
both systems are in the same universality class. To settle this question numerically, an algorithm similar to that for 
uncorrelated percolation is needed to generate a random vortex tangle, consisting of intertwined closed vortex loops. 
Such an algorithm is provkled by the Vachaspati-Vilenkin algorithmc3 developed to study defect formation in rapid 
(cosmic) phase transitionsEj. r. 

Quantized vortex lines emerge when a U(l) symmetry is spontaneously brokenu. The vacuum of a state with such 
a broken symmetry is characterized by a phase factor exp(z<y9), with the phase tp changing by 2tt when a vortex of 
unit topological charge is circled once. 

Consider discretizing both space (to a three-dimensional square lattice, say, with N sites) and also the vacuum 
manifold, so that if takes only certain discrete values in the interval < ip < 2tt. The restriction to this interval 
assures that only vortices of unit strength are generated. Imagine visiting a site of the lattice and randomly assigning 
it one of the discrete phase values. The random choice guarantees that phases at different lattice sites are uncorrelated. 
The Vachaspati-Vilenkin algorithm, with ip taking the minimum number three of different values, can be summarized 
in pseudocode as: 

while i < N 



vortex 




• • 

2 4 ^* 

4 



FIG. 6: Illustration of Kibble's geodesic ruIeO. Top: Plaquette pierced by a vortex. The path segments labeled by "a" show the 
path taken by the phase along the vacuum manifold (circle at the right) when one moves from the left-lower to the right-lower 
corner on the plaquette circling the vortex. Bottom: Although all three available values of the phase appear, no vortex pierces 
the plaquette. When the right-lower corner is reached, where the phase is 47r/3, the phase moves back to 2n/3 rather than to 
to complete the circle, representing the vacuum manifold. 
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FIG. 7: Illustration of two vortices entering a (three-dimensional) unit cell. 



generate ifi S 
i + + 
end while 



{0, 



2^ 



After all lattice sites have been visited in this way, vortices are traced by going around each plaquette of the lattice 
once. When upon returning to the starting site, a phase difference of 27r is found, it is concluded that a vortex 
penetrates the plaquette. |-in going from one lattice site with the phase ipi to a neighboring site with the phase 
if2. Kibble's geodesic ruleEZi is used to select the [Shortest path along the circle representing the vacuum manifold 
connecting (pi and ip2 (see Fig. ^ for two examples^). 

When two vortices are found to enter a unit cell (see Fig. |^), the two incoming and two outgoing vortex segments are 
randomly connected. After all the plaquettes of the lattice have been considered, a random vortex tangle, consisting 
of closed loops can be traced out. p. 

The fractal dimension of the vortex tangle found iiLthe initial simulationsE^ were consistent with that of Brownian 
random walks, D = 2. This result can be understoodD by realizing that a dense random vortex tangle was simulated, 
where a single vortex experiences an effective repulsion from the neighboring vortex segments as the volume occupied 
by these segments is no longer available. This is similar to a polymer in a dense solution, which also executes a 
Brownian random walk. 

In the Vachaspati-Vilenkin algorithm, the allowed phase values have equal probability so that the vortex tangle 
generated has always the same average vortex line density. To be able to vary the density, VachaspatiB lifted the 
degeneracy of the vacuum manifold by giving one phase value a bias, P27r/3 — Piir/s 7^ Pa- This can be achieved by 
choosing for examplaJ 



P27T/3 — P4tt/3 



^Ce"/^, Po = Ce- 



(36) 



with = 2exp(77/2) -|-exp(— 77). Since all phase values are needed for a vortex to pierce a given plaquette, the bias 
obstructs the formation of vortices and decreases the vortex line density. The bias parameter rj plays a role similar to 
that of the occupation probability p, determining the cluster density in percolation theory. 

Imagine starting in the regime with only small vortex loops present and gradually changing the bias so that more 
and at the same time larger loops appear. It turns out that for small enough values of rj an "infinite" vortex loop is 
frequently generated during updateaJ. The division between finite and infinite vortex loops is,, of course, somewhat 
arbitrary on a finite lattice. A possible choice is to classify vortices with a length larger thariEa 



d 

— -Z 

It: a 



(37) 



as infinite, where L is the lattice size. The length scale Lc is the average length of Brownian random walks traversing 
the lattice once. In other words, most infinite loops wrap themselves around the lattice and thus have a nonzero 
winding number. Universal quantities such as the critical exponents should be independent of the precise choice. On 
an infinite lattice, a threshold r] = rjc exists, above which infinite vortex loops are absent and large, but finite loops 
exponentially suppressed. The sudden appearance of an infinite vortex Ipiop as the threshold 77 = 77c is approached 
from above is similar to that of a percolating cluster in percolation theorytl. For an infinite vortex to appear, the line 
tension 6 needs to vanish, so that the Boltzmann factor in the loop size distribution (35) becomes unity and no longer 
suppresses large loops, i.e.. 



(38) 
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In analogy with percolation theory, cf. Eq. ([l|), the vanishing of the line tension is specified by an exponent a, 

cx (ry - 77e)i/^ (39) 

At the threshold, vortex loops proliferate in the same way that clusters do at the percolation threshold. 

To characterize the transition, the same definition of critical exponents can be used as in percolation theory, leading 
to the same scaling relations ( p^ found there. Numerical simulations of the transition that random vortex tangles 
undergo when changing the bias using this algorithm showed that the values for the critical exponents are more or 
less consistent with those of three-dimensional uncorrelated site percolationE3: 



critical 


percolation 


minimally 


exponent 


theory 


discretized U(l) 


a 


0.45 


0.46(2) 




0.41 


0.54(10) 


7pcr 


1.80 


1.59(10) 



Other values obtained in simulations with the discretized vacuum manifold given more than the minimal three point: 




show even better agreement with percolation theory, implying that uncorrelated percolation and random vortex tangles 
are indeed in the same universality class. In other words, the close similarity in describing both phenomena discussed 
above extends to the level where fluctuations, determining the numerical values of the critical exponents, are included. 



C. Correlated Vortex Tangles 



In Sec. [IB, the extension of the theory of uncorrelated percolation to describe thermal phase transitions in corre- 
lated systems was discussed. It was concluded that correlated and uncorrelated percolation differ only in the specific 
values of the exponents a and r, specifying the cluster size distribution (^. The general framework applies to both 
cases, with the various expressions giving the critical exponents in terms of a and r remaining unchanged as they are 
protected by scaling laws. 

Here, a similar extension is discussed for the random vortex tangle to describe superfluid phase transitions. An 
arbitrary — not just random — vortex tangle is again described by the loop size distribution (p5|)cj. In the superfluid 
phase, large vortex loops are exponentially suppressed. When the temperature increases, larger loops appear and the 
vortex tangle becomes denser. As the critical temperature is approached, the line tension 9 vanishes, thereby allowing 
vortex loops to proliferate and thus to disorder the superfluid phase. The vanishing of the line tension is specified by 
the exponent a, 

e (X (Tc -r)^/'". (40) 

As in correlated percolation, cr can be related to the fractal dimension D, introduced in (p9|), leading to the same 
relation (H) . The value of the fractal dimension differs in general from the value obtained for a random vortex tangle 
or for Brownian random walks {D = 2). At criticality, the correlation function ( |3^ ) shows algebraic behavior, 

G(-) - (41) 

as does the correlation function (^ defined in percolation theory. Using scaling laws, one can again express all the 
critical exponents in terms of the exponents a and r this time specifying the loop distribution. 

However, with the exception of random vortex tangles discussed in the preceding subsection, usually not the 
percolation set of critical exponents, but a different set is used to describe vortex loops. The difference can be seen 
by considering the correlation function in percolation theory [above Eq. (^)] and for random walks [Eq. (|3^)], or more 
precisely their sum rules: 

Gporix) cc (X Y s'"" (42) 

X s s 

and 

Y G{x) cx Y « e-''^"", (43) 
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FIG. 8: CoHiputed vortex tangle in pure superfluid turbulence without normal component. From Araki, Tsubota, and 
NemirovskiitJ. 



respectively, where the sum is over all lattice sites. The summands at the right sides of Eqs. (^) and (^) are 
seen to depend differently on the size (s and n, respectively). To understand the origin of this difference, note that the 
left sides equal the susceptibility x for the respective order parameters. In general, x is calculated from the partition 
function by differentiating it twice with respect to the conjugate field coupling linearly to the order parameter. The 
conjugate field is, for example, the magnetic field for magnetic systems, and the pressure for liquid-gas transitions. 
In other words, these correlation functions are tied to the definition of the order parameter, each having its own set 
of critical exponents. 

For random walks, the set of critical exponents read in terms of the exponents a and r specifying the loop size 
distributionme: 



2cr 



7=-, 11 = '^- D. 
a 



(44) 



Here, the critical exponent which specifies the vanishing of the order parameter when the critical temperature is 
approached, is given a subscript "ce" (for critical exponent) to distinguish it from the inverse temperature. These 
formulas should be compared to the analogous formulas in En. (|lO|) for percolation. The expressions for 7per and 7 
are directly related to the size-dependences in Eqs. ( ^2|) and (^3|). Brownian random walks, having fractal dimension 
D — 2, separate self-avoiding vortex loops {D < 2) from self-seeking ones {D > 2). 

The two different sets of critical exponents, for percolation and for random walks, both satisfy the scaling laws. As 
only one correlation length features in both descriptions, cf. Ref. ^3|, the two sets have the exponent v in common, 
implying that they are not independent. Consider, for example, the 0(2), or XY-model, whose thermal phase 



transition admits a description in terms o f bot h correlated percolation and vortex loops. Let (t„ 
cluster size distribution discussed in Sec. [IB, and cr and r the vortex loop size distribution (p. 



and Tpor specify the 
, then 



T — 1 (T 
'^pcr 1 ^per 



(45) 



or written in terms of fractal dimensions. 



D a„ 



(46) 



where Dpor denotes the fractal dimension of the (bond) clusters, while D denotes that of the vortex loops at criticality. 



D. Generating Vortex Tangles 

The vortex tangles in superfluids discussed so far were generated by thermal fluctuations in the vicinity of the phase 
transition. At least two other ways to generate vortex tangles in superfluids are known. 

(i) Consider setting up a counterflow, so that the superfluid and normal components move with different velocities, 
or dragging an object such as a grid through the superfluid. When a critical velocity is reached, Xhe superfluid 
becomes turbulent through the formation of a random vortex tangle made jHC of quantized vortex loopso (see Fig. |^) . 
For a constant, uniform counterflow, numerical simulations by SchwarzEJ showed that the vortex tangle can be 
self-sustained. In such counterflow turbulence, the normal component is possibly laminar and only the superfluid 
component turbulent. In superfluid grid turbulence, on the other hand, the superfluid and normal components are 
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believed to become locked together and behave like in classical turbulence, at least on scales where the quantization 
of vortices is irrelevant. 

(ii) Consider heating a small region of the superfLuid by, for example, radiating the cell with neutrons as in the 
Helsinki experiment on superfluid ^He-B in rotatioiO. When a neutron is absorbed by the supcrfluid, the energy 
deposited causes a small region to briefly heat up to the normal phase. After the subsequent rapid cool-down, jO, vortex 
tangle emerges. This scenario for defect formation after a rapid phase transition was first proposed by KibbletZI in the 
context of cosmology and extended to condensed matter by Zurekcj. From the present perspective, the emergence of 
a vortex tangle after a rapid temperature quench is not so much the formation of defects, as well as the "freezing" into 
permanency of the dense vortex tangle rendering the hot bubble normal. In this sense, the vortex tangle emerging 
after the quench is a remnant of the normal phase. In addition to infinite loops, the normal phase has also many small 
vortex loops, which initially survive the rapid temperature quench. But as time passes, these small loops shrink to 
zero, leaving behind only larger ones traversing a macroscopic portion of the system (see Fig. The larger loops by 
no means remain unchanged in time as they, through reconnections, spin off small loops, or hook up with other vortex 
loops to grow. The final configuration depends, of course, on the conditions under which the quench is performed, 
such as the presence of superflow. Apart from infiating vortex loops, which explains the expansion of the vortex tangle 
in Fig. ^, a superflow can also lead to additional mechanisms £ax defect formation, vortex loops wrapping around the 
bubble being created to screen the superflow from the interiorcZI. 



IV. VORTEX ACTION 



The discussion of vortex loops so far was qualitative as few details are needed to describe phase transitions in terms 
of vortex proliferation. In this section, a more detailed description of vortex lines is given. 



A. Schrodinger Equation 

To obtain the action describing vortex lines in a superfluid, the differential equation ( p7[ ) satisfied by the transition 
probability Ar„(x x') is written in the suggestive form 

dsKsix ^ x') = iftV2if,(x ^ x') (47) 

by introducing the arclength parameter s. 

In this way, the differential equation becomes analogous to the time-dependent Schrodinger equation for the "time" 
evolution operator iirs(x — > x'), satisfying 

V'(x',s)- y"d''xX,(x^x')V'(x,0), (49) 

where '0(x, s) is the wave function at positi(m,x and "time" s. The arclength parameter does not represent real time, 
but is the Schwinger propertime parameterEB. It parametrizes the path taken by the random walker, which is now 




FIG. 9: Computed evolution of a vortex tangle in superfluid ^He-B after a rapid temperature quench of a heated bubble in 
the presence of a superflow (perpendicular to the final configuration, resembling a snow chain). From Aranson, Kopnin, and 
VinokuiO. 
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FIG. 10: Interaction between two line segments belonging to two different vortex loops. 



pictured as a (bosonic) particle trajectory. The corresponding action reads: 



(50) 



where an extension to arbitrary many particles (labeled by q) is given, and the second term resulting from the 
modification (|3^) is included, with 9 the vortex line tension. The action is given a subscript to indicate that 
the interaction between vortex lines has been ignored for now. With one of the spatial coordinates interpreted as 
the time coordinate, tte action (^^ describes for 6* > the worldlines of relativistic particles of mass •\/2i? in d- 
dimensional spacetimeEll. The particles are at zero temperature, with the quantum fluctuations representing the 
thermal fluctuations of the vortex system at finite temperature. 

To transcribe the (Euclidean) action 
the dimcnsionless combination S/h is to 
Boltzmann factor: 



, obtained using the quantum-mechanical analog into an Hamiltonian, 
DC replaced with fiH , so that the exponential factor exp(— 5/fi.) becomes a 



exp(-S'/n) exp(-/3ff). 



(51) 



The Planck constant h characterizing the quantum fluctuations of the equivalent zero-temperature particle system is 
replaced with the inverse temperature /3 characterizing the thermal fluctuations of the vortex system. The transcription 
is achieved by the change of variable 



hhf3s, 



where the factor i is included for mere convenience. The transcribed action then reads 



2de 



(52) 



(53) 



where the prime on the arclength parameter is suppressed, and the parameter cl, which is determined by the vortex 
line tension 0, has the dimension of an energy per unit length. It follows from this expression that the line tension 
should scale linearly with a to obtain a flnite value for the parameter cl in the continuum limit. The correlation 
function (pa) reads in these variables 



Go(|x-x'|) = ^ / dsX,(x^x') 



^^rd,sm 



d/2 



\47rs / 



exp 



/3(x-x')^ 
4 ~s 



(54) 



The exponent d/2 in the first factor of the integrand at the right hand expresses the fact that the random walker can 
roam in d dimensions, each giving a contribution yf(jlAi^. 



B. Interaction 



1 . Superfluids 



In this subsection, the interaction between vortex loops in a superfluid is considered. The pairwise interaction is 
mediated by the Goldstone mode associated with the spontaneously broken U(l) symmetry. Being gapless, this mode 
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FIG. 11: Suppression of vortex loops and creation of open vortex lines when an external magnetic field is applied. 



leads to a 1/R interaction in three space dimensions, 

(27r)3 A:2 ^ 



I 



(55) 



Each segment of a vortex line interacts via this potential with all other line segments of the same vortex as well as of 
the other vortices (see Fig. |l^). The potential energy y of a vortex tangle is therefore of the form 

^ = ^f^E/ds/d^'^«(^)-^^?'(^')- (56) 

9,9' 

Here, g is the charge of a vortex given by 

(57) 

with h/m the circulation quantum and ps the superfluid mass density, while R is the distance between two line 
segments, R = |xg(s) — Xq'(s')|. The mass parameter m denotes the mass of the atoms (for bosonic superfluids) or 
that of a pair of atoms (for fermionic superfluids). The integrations in Eq. ( p6| ) arc along the vortex loops parametrized 
by their arclength parameters s and s'. Since V decreases with increasing i?, the pair interaction between vortex 
segments with like-charges is repulsive. The total Hamiltonian H describing a vortex tangle consists of the sum of 
the free part ( p3| ) and the potential energy (|5^), 

H = Ha + V. (58) 

For two slightly deformed rectilinear vortices parallel to the z axis, each of length L and separated by a distance 
r{z), the right hand of Eq. ( |5^ ) reduces to the well-known logarithmic form 

V='^l\zl" d.- , ' = In . (59) 



471- Jo Jo ^r^{z) J- -^0 



The dependence on the vortex length L in the argument of the logarithm derives from the slow fall-off of the 1 /R 
interaction. 



2. Superconductors 

The results obtained for superfluids are easily adapted for superconductors. The major difference is that the 
would-be Goldstone mode turns the long-ranged electromagnetic interaction into a screened one through the Higgs 
mechanism. This leads to the modification 
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dz' \ 



dz 



FIG. 12: Two line segments of the same vortex (solid line) interacting via the 1/ R potential (dashed line). 



in Eq. (p6|), where A is the London penetration depth. The charge g of a magnetic vortex is again given by Eq. (|57] ) 
for a superfluid, although it is usually put in the equivalent form 

9=—' (61) 

where $0 is the magnetic flux quantum. Note that g is independent of the electric charge e. The e-dependences in 
<i>o and the penetration depth precisely cancel in the right hand of Eq. (61). Since the interaction is screened, the 
potential energy between two parallel slightly deformed rectilinear magnetic vortices separated by a distance r(z) is 



with Kq a modified Bessel function of the second kind, which for small argument x behaves as Kq{x) ^ — ln(x). The 
role of infrared cutoff is now played by the London penetration depth A. 



C. External Field 



From the present perspective, the major effect of turning on an external magnetic field is to suppress fluctuations 
along the field. In other words, the random walker representing the vortex trajectory is forced to move only in 
the direction of the applied field, and cannot move backwards (see Fig. ^l]). The external field thus suppresses the 
formation of vortex loops and (under certain conditions discussed below) instead leads to open vortices along the field 
direction, starting at the bottom of the sample and terminating at the top. The vortices can now only fluctuate around 
their equilibrium position by making excursions in the plane perpendicular to the field. In the quantum-mechanical 
analog, this corresponds to taking the nonrelativistic limit. 

To study the effect of applying an external field, which is assumed to be parallel to the z axis, write 

(x-x')' =r2+c2z^ (63) 

where r(z) = [x{z),y{zy\ are the coordinates perpendicular to the field, parametrized by z. The parameter c is 
introduced (via the replacement z — > cz) to facilitate interpreting the spatial coordinate z as a, time coordinate. The 
observation that fluctuations along the field are suppressed is implemented by letting c — > cxd. In this-Jimit, the 
integral over the Schwinger propertime parameter s in Eq. (|5^) can be approximated by the saddle pointEHl 

.^^z, (64) 
connecting s to the timelike variable z, with the result in three dimensions 

^Go(lx-x'l) « (A^'' e-/5e.(.+rV2.). (65) 

Here, the parameter c is scaled away again by the replacement cz ^ z thereby undoing the earlier replacement in Eq. 
(|63|). The exponent in the prcfactor at the right hand is unity now instead of 3/2 as it was in the absence of a field. 
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It indicates that after applying the magnetic field only the t wo dire ctions perpendicular to the field are still available 
to the random walker to roam, each giving a contribution I3/Aitz. The Hamiltonian ( ^3|) becomes in this limit 

ffo - 5] / ds [i x2(s) +el]^Y.h^ [5 ^Lr^(^) + el] , (66) 

where use-Js made of Eq. (|6j). This procedure of taking the nonrelativistic limit is equivalent to the one used 
by FetterLj of expanding in small displacements in the xy plane. With z interpreted as the time coordinate, this 
Hamiltonian becomes, using the transcription (pT]), an action describing nonrelativistic bosons of mass eL moving in 
a plane in the presence of a constant background potential eL- As in the absence of an external field, the vortices, 
which are now open lines, are interpreted as the worldlines of particles, with the quantum fluctuations representing 
the thermal fluctuations of the vortex system. The limit taken in Eq. ( |66| ) precisely corresponds to the nonrelativistic 



limit of the relativistic Hamiltonian (53). The interaction between two "nonrelativistic" vortex lines is given by Eq. 
( |6^ ) as the case of slightly deformed rectilinear vortices considered there precisely describes vortices in the presence 
of an external field. 

Not only two line segments belonging to different vortices interact via the 1 /R potential, but also two line segments 
of the same vortex (see Fig. Ol), giving the contribution to the energy 



^ dz dz' ^ = I- / dz In ^ , (67) 

4^ Jo Jo ^J[z-z'Y 4^io \i) 

where the coherence length ^ is taken as ultraviolet cutoff. This contribution amounts to a renormalization of the 
parameter el- In the following, that parameter, introduced in Eq. (^3|), is assumed to be given by the integrand at 



the right hand of Eq. (67), representing the self- interaction energy per unit length, i.e., 

eL-l^ln(Ac), (68) 

where k — A/^ is the Ginzburg-Landau parameter. The total Hamiltonian iJ describing vortices in the presence of 
an applied magnetic field is then given by the sum of the free part (|66|) and the interaction part 

-Hint — -T' 



^1 i,Ko (^i^iM^iVM^ _ci,„^ J dz/iext[r,(z)], (69) 



where the last term represents the interaction of the vortex lines with the external field. This field penetrates the 
sample only if the sum of the constant terms in Eqs. (|6|) and ( |69| ) becomes negative, i.e., when 

hcxt > = ^ci , (70) 
<Po 

defining the lower critical field /id- Below this critical field, the system is in the Meissner state where the external 
field is excluded from the superconductor. Above hc^ , the system is in the mixed state where the external field invades 
the superconductor through an array of magnetic vortex lines. When fluctuations CEwtbe ignored, the repulsive pair 
interaction between the vortices makes them form a triangular Abrikosov flux latticeEJ. 



V. ENTANGLED VORTEX LINES 



In this section, vortex tangles of a different kind than the ones above are discussed. Instead of consisting of closed 
vortex loops, the tangles considered here consist of open vortex lines, connecting the top and bottom of the sample 
due to an applied magnetic field (for superconductors) or rotation (for superfluids). 



A. Vortex Lattice Melting 



The equivalent particle system of nonrelativistic interacting bosons in 2+1 dimensions (with the z coordinate 
interpreted as aptimelike variable) specified by the Hamiltonian H = Ho + i^int given in Eqs. ( |6^ ) and ( |69| ) was 
taken by NelsonEa as starting point to investigate the effect of thermal fluctuations on the Abrikosov vortex lattice. 
A similar approach to study fluctuations in vortex lattices in rotating ^He and in superconductors was pioneered 
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FIG. 13: Artist's impression of the Abrikosov flux lattice (left panel) and the entangled vortex liquid (right panel). Periodic 
boundary conditions in the z direction are assumed. 

by FetterEi. In the particle system, the vortex lattice corresponds to a Wigner crystal resulting from the repulsive 
pairwise interaction between the particles. The thermal fluctuations of the vortex system translate into quantum 
fluctuations of the equivalent particle system, which, when violent enough, were predicted to melt the crystal. The 
resulting phase in the equivalent particle system is the supcrfluid phaseE3, which is generally believed to be the sole 
alternative to the crystal phase for interacting bosons at low temperatures in a clean system without impurities. The 
superfluid phase corresponds to entangled vortex lines in the vortex system (see Fig. p^ ). 

The importance of quantum fluctuations in a system of interacting bosons is measured by the so-called de Boer 
parameter AE3. This dimensionless parameter can be obtained by writing the kinetic term at the right hand of Eq. 
( |66| ) and the pair interaction term in Eq. ( |69|) as 

(71) 

where u = fiag^zjl-n and 

= /3aVg2eL/27r, (72) 

with all length scales now being given in units of a typical length a, e.g., the lattice spacing when working on a 
lattice. The vortex lattice melting transition into an entangled vortex liquid was nipHerically shown to take place at 
the critical value Am ~ 0.062E3. There is ample numericalE^I as well as experimentelEj evidence that the transition is 
discontinuous in a clean system. Further experimental results have been reportedEj, indicating that the discontinuous 
melting line in the h-T phase diagram terminates at a critical point in the high-field, low-temperature corner. 
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B. Feynman's cooperative exchange rings 

As was pointed out by Nelsonll, the entangled vortex liquid can be elegantly undeijstood in the equivalent particle 
system, using Feynman's cooperative exchange ring theory of the A transition in "^Ha^j. Bose-Einstein condensation 
(BEC) is understood in this picture as follows. At finite temperature, the (3-|-l)-dimensional ^He particle trajectories 
form closed loops because the (Euclidean) time coordinate is compactified, taking values only in the interval [0,hp]. 
For an ensemble of bosons, the boundary conditions are periodic, meaning that the configuration at time and at time 
hp are identical. In the normal phase, most bosons execute during this time interval a random walk in space (see left 
panel of Fig. |lj), which starts and ends at the same position (indicated by a dot). It means that the particles, being 
distinguishable, behave classically. When the critical temperature is approached from above, the particles lose their 
identity and become indistinguishable. In Feynman's theory this is reflected-by the formation of so-called cooperative 
exchange rings, where individual worldlines hook up to form larger loopsEj (see right panel of Fig. |lj). A particle 
in such a composite ring moves in imaginary time along a trajectory that does not end at its own starting position, 
but at that of another particle. Hence, although the initial and flnal configurations are identical, the particles in a 
composite ring are cyclically permuted. 

To determine the condition for this to happen, consider an ideal Bose gas for simplicity. The non-interacting 
particles are described by a time evolution operator ii'r(x — > x'), satisfying the time-dependent Schrodinger equation 
( ^ ) with the Schwingcr propcrtime parameter s introduced in Eq. ( |4q ) now replaced by the (Euclidean) time 



(73) 



20 



y,z 

A y,z 




FIG. 14: Trajectories of non-mteracting bosons in three space dimensions executed during the (imaginary) time interval h/S. For 
ease of representation the y and z axes are lumped together. The starting points of the Brownian random walks are indicated 
by dots. In the normal phase (left panel), most random walks end at their own starting point, whereas in the Bose-Einstein 
condensed phase (right panel) large cooperative exchange rings are formed, with the trajectory of one particle ending at the 
start ing point of another. In this way, particles in a ring are cyclically permuted after an imaginary time ft/3. (Adapted from 
Ref. |6|.) 



The right hand explicitly depends on the particle mass m. It follows that during the time interval At — fi-/3, the 
particles execute a Brownian random walk with an average length 

Ln,-^'^ (74) 

determined by the de Broglie thermal wavelength Xp. If the length scale Lnp becomes, upon reducing the temperature, 
on the order of the interparticle distance, individual worldlines can hook up to form large loops. The particles 
participating in the exchange rings become indistinguishable and the system Bose-Einstein condenses. When worked 
out, this condition gives an estimate of the critical temperature T^. 

A ring made up of the worldlines of w particles wraps around the imaginary time cylinder w times. BEC correspoiids 
in this picture to the appearance of long loops, wrapping arbitrarily many times around the imaginary time cylinders. 



In other words, as for vortex loops before (see Sec. |III|), this transition corresponds to a proliferation of loops — not 
of vortex loops, but of worldline loops this time. BEC can therefore be described as a percolation processEj. The 
worldline loop size distribution^^ 

C oc w-"e-("-"")'", (75) 

with a = ^Pfi and t = d/D + 1, gives the density of rings containing w particles, or, equivalently, the density of loops 
wrapping around the imaginary time cylinder w times. The "line tension" a — ac of the worldlines is determined by 
the chemical potential /i [as was the Boltzmann factor in the cluster size distribution (^4|) in Mayer's theory of the 
liquid-gas transition] and vanishes when the critical temperature is approached from above [cf. Eq. (|l8|)]. Specifically, 

a-ttc oc (T-Te)'/". (76) 

The worldline loop size distribution then changes from an exponential decay with increasing loop length to an algebraic 
decay. The winding number w in Eq. ([75|), denoting the number of particles contained in a ring plays the role of the 
number n of steps in the loop size distribution (p5[). In other words, a straight path segment traversed during a single 
step by a classical random walker now corresponds to a random trajectory traversed by a quantum particle during the 
time interval h(3 [with an average length Lnp given by Eq. (|^) if the particle is non-interacting] . Since the loop size 
distribution is of statistical nature, details which req uire a resolution greater than the step size are irrelevant — whence 
the similarity in the loop distributions ( [ts]) and (^). Given the two exponents a and r specifying the worldline loop 
size |distribution ([75[), the critical exponents characterizing the phase transition can again be obtained using scaling 
lawso 

Feynman's picture of the A transition is dual to the one based on vortex proliferation discussed earlier. Depending 
on whether the critical temperature is approached from below or above, vortices (in the superfluid phase), or worldlines 
(in t he n ormal phase) proliferate. This is similar to the duality discussed in the context of percolation at the end of 



Sec. [I A, where clusters of the minority phase are considered in the sea made up of the majority phase: here either 



vortices in a sea of proliferated worldlines, or worldlines in a sea of proliferated vortices are considered. 
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Translated back to the vortex system, Fig. |l^ represents the projections of the vortex hnes onto a plane perpendicular 
to the applied field. The periodic boundary conditions of the equivalent particle system are somewhat unnatural for 
the vortex system, but are believed not to dramatically change the results, provided L, playing the role of h(3 in 
the equivalent particle system, is large enough. The left panel of Fig. |lj mimics the Abrikosov flux lattice (apart 
from the lattice symmetry) , with the vortex lines making only small excursions from the straight lines connecting the 
endpoints, while the right panel mimics the entangled vortex liquid. Precisely such snapshots of vortex positions were 
recently obtained in numerical simulations by Sen, Trivedi, and CeperleyE3. 



C. Order Parameter of Entangled Vortex Liquid 



A superfluid phase is characterized by the superfluid mass density ps giving the response of the system to an 
externally imposed velocity v. A computationally simple, yet powerful way of obtaining the|-superfluid mass density 
on a finite lattice with periodic boundary conditions was introduced by Pollock and CeperleyE3, who related ps to the 
winding of particle trajectories around the lattice (not in the time direction, but in the space directions). 

The free energy F = —f3~^ ^^{Z) acquires an additional term after the boost 

A^^„ = iVpsV^ (77) 
In the equivalent particle picture, this boost is implemented by letting 

dr 

imr^ \m[rq - v)^ w imr^ - "^"d"^ ' 

leading to 

e-^^^" = (exp(i^w.v)), (79) 

with the winding vector: 

9 







To understand the appearance of the imaginary unit at the right hand of Eq. (|7^), note that in real time, particles 
are described by the phase factor exp(iS'/?i). When going over to the Euclidean time variable r, the term linear in 
the velocity retains the imaginary unit in the phase factor. 

Since the integrand in the definition of the winding vector is a total derivative, only the endpoints of the integration 
contribute. In numerical simulations on a finite lattice, usually also periodic boundary conditions in the spatial 
directions are imposed. The winding vector then becomes 

w-^[r,(n/3)-r,(0)] + Ln, (81) 

9 

with L the linear size of space and n a d-component vector with integer components, denoting the number of times 
the trajectories wrap around space in the x, y, ■ ■ ■ direction, respectively. Because of the periodic boundary condition 



in the r direction, the sum in Eq. (81) is always zero, even when particles are cyclically permuted, so that a nonzero 
value for the winding vector obtains only when worldlines wrap around space, represented by the last term in in Eq. 
(pl|). In the language used in connection with the length scale (|37|), such loops axe classified as infinite. Expanding 



both sides of Eq. ( |79D to quadratic order in the applied velocity v, one finds thatt 

(w2) = dNPg\ (82) 



with g the vortex charge introduced in Eq. (pT)), = {h/m)'^ps. It shows that the superfluid mass density is nonzero 
when infinite worldline loops are present. That is, Ps directly signals the proliferation of these loops and is therefore 
the proper order parameter describing the entangled vottex liquid. 

The numerical simulations by Nordborg and BlatteitJ show that the superfluid mass density sharply rises at the 
critical value Am ~ 0.062 from ps ~ in the crystal phase to ps ~ p in the superfluid phase, where p is the total mass 
density. 
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VI. CONCLUSION 



In these lecture notes, thermal phase transitions taking place in superfluids were discussed from a geometrical 
perspective. Two seemingly different geometrical approache&,were treated: one modeled after cluster percolation and 
the other after loop proliferation first proposed by Onsageiu. Apart from peculiarities resulting from differences in 
definition, the two geometrical approaches were, however, seen to be two faces of the same underlying formalism. 
Central in the description was the size distribution of clusters and loops, specified by two exponents from which in 
turn all the critical exponents characterizing the superfluid phase transition follow through scaling relations. At the 
critical point, two different types of clusters percolate and also two different types of loops proliferate, depending 
from which side the transition temperature is approached. The percolating cluster always belongs to the minority 
phase which then turns into the majority phase. The two percolation descriptions are dual to each other in that both 
phases are interchanged. Also the two loop descriptions are dual to each other, with finite vortex loops featuring 
in the superfluid phase and finite worldline loops featuring in the normal phase. These finite loops proliferate when 
approaching the critical temperature after which their roles are interchanged. 

Finally, the similarity between the cluster distribution (|8|) and the worldline loop distribution (|7|), both depending 
on the fugacity exp(— a), indicates that the descriptions in terms of clusters and worldline loops are in fact closely 
related. This connection is reinforced by the surprising presence of the de Broglie wavelength (which depends on h) 
in the classical partition function (13). 
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